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ABSTRACT 

Several  publications  describe  research  carried  out  at 
NLR  on  the  thermal-gravitational  modelling  and  scaling 
of  two-phase  heat  transport  systems  for  spacecraft 
applications.  They  dealt  with  mechanically  and  capillary 
pumped  two-phase  loops.  The  activities  pertained  to 
pure  geometric,  pure  fluid  to  fluid,  or  hj^rid  scaling 
between  a  prototype  system  and  a  model  at  the  same 
gravity  level,  and  between  a  prototype  in  micro-gravity 
and  a  model  on  earth.  Recent  publications  also  include 
the  scaling  aspects  of  a  prototype  loop  for  a  Moon  or 
Mars  base  application  and  a  terrestrial  model. 

The  work  discussed  here  was  carried  out  in  the  last 
couple  of  years.  It  concerns  scaling  to  super-gravity 
levels,  and  was  done  because  a  promising  super-gravity 
application  for  (two-phase)  heat  transport  systems  can 
be  the  cooling  of  high  power  electronics  in  spinning 
satellites  and  in  military  aircraft.  In  such  aircraft,  power 
electronics  can  be  exposed,  during  manoeuvres,  to 
transient  accelerations  up  to  say  12  g  (120  m/s^).  The 
discussions  include  pulsating  two-phase  ioops  and  also 
oscillating  (pulsating)  heat  transfer  devices,  for  which 
several  important  applications  in  different  acceleration 
environments  have  been  identified. 

INTRODUCTION  AND  SCALING  APPROACH 

Many  examples  of  the  scaling  of  two-phase  flow  and  heat 
transfer  can  be  found  in  the  power  industry  and  in  the 
process  industry.  The  scaling  of  physical  dimensions  is  of 
principal  interest  in  the  process  industry,  where  large-scale 
industrial  systems  are  studied  by  reduced  scale  laboratory 
model  systems.  The  scaling  of  the  working  fluid  is  of 
principal  interest  in  the  power  industry.  There,  large-scale 
industrial  systems  (characterised  by  high  heat  fluxes, 
temperatures,  and  pressures)  are  translated  in  full  size 
systems,  operating  at  lower  temperature,  heat  flux  and 
pressure  levels  (e.g.  the  scaling  of  a  high  pressure  water- 
steam  system  by  a  low  pressure  refrigerant  system  of  the 
same  geometry). 


The  thermal/gravitational  modelling  and  scaling  of  two- 
phase  heat  transport  systems  was  carried  out  at  NLR: 

-  For  a  better  understanding  of  two-phase  flow  and  heat 
transfer  phenomena. 

-  To  provide  means  for  comparison  and  generalisation 
of  data. 

-To  develop  a  useful  tool  to  design  two-phase  systems 
and  components,  in  order  save  money  and  to  reduce 
costs. 

The  main  objective  of  scaling  space-oriented  two-phase 
heat  transport  systems  and  system  components  is  the 
development  of  reliable  spacecraft  systems,  of  which  the 
reduced  gravity  performance  can  be  correctly  predicted 
using  results  of  terrestrial  experiments  with  scale  models. 

Scaling  of  spacecraft  systems  can  be  useful  also: 

-  For  in-orbit  technology  demonstration.  The  performance 
of  spacecraft  heat  transport  systems  (prototypes)  can 
be  predicted  using  the  outcomes  of  in-orbit  experiments 
on  models  with  reduced  geometry  or  different  working 
fluid. 

-To  define  in-orbit  experiments  to  isolate  a  typical 
phenomenon  to  be  investigated  (e.g.  the  exclusion  of 
gravity-induced  disturbing  buoyancy  effects  on  alloy 
melting,  diffusion  and  crystal  growth)  for  a  better 
understanding  of  the  physical  phenomena. 

The  magnitude  of  the  gravitational  scaling  varies  with  the 
objectives,  from: 

-  1  g  to  1 0’®  g  (random  direction)  for  the  terrestrial  scaling 
of  orbiting  spacecraft. 

-  1  g  to  0.16  g  for  Moon  base  systems,  to  0.4  g  for  Mars 
base  systems. 

-  10’^  g  or  lO"®  g  to  1  g  for  isolating  gravity-induced 
disturbances  on  physical  phenomena  investigated. 

-  Low-g  to  another  or  the  same  low-g  level  in  low-g 
aircraft  or  sounding  rockets. 

-  1  g  to  say  12  g,  for  electronics  cooling  in  military 
aircraft,  characterised  by  (transient)  super-gravity 
conditions. 
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NLR  developments  originally  pertained  to  the  scaling  of 
mechanically  and  capillary  pumped  two-phase  loops  for 
use  in  microgravity  (Refs.  1  to  7).  The  activities  were 
extended  to  applications  in  Moon  and  Mars  bases  (Ref. 
8).  Modelling  and  scaling  activities  of  the  last  two  years 
focus  on  mechanically  and  capillary  pumped  two-phase 
loops  and  oscillating  heat  transport  devices  for  use  in 
transient  super-gravity  environments,  encountered  during 
military  aircraft  manoeuvres.  Similar  activities  are  reported 
to  be  currently  carried  out  elsewhere  (Refs.  9, 10). 

It  is  recalled  from  the  aforementioned  publications  that  even 
in  single-phase  systems,  scaling  is  anything  but  simple, 
since  flow  and  heat  transfer  is  equivalent  in  the  model  and 
the  original  (prototype)  only  if  the  corresponding  velocity, 
temperature  fields  and  pressure  fields  are  identical. 
Dimensionless  numbers  can  be  derived  either  from  the 
consen/ation  equations  (mass,  momentum  and  energy)  or 
from  similarity  considerations,  based  on  dimension 
analysis.  The  aforementioned  identity  of  velocity, 
temperature  and  pressure  fields  is  obtained  if  all 
dimensionless  numbers  are  equal  in  model  and  prototype. 

Scaling  two-phase  systems  is  far  more  complicated  as: 

-  In  addition  to  the  above  fields,  the  spatial  density 
distribution  (void  and  flow  pattern)  is  to  be  considered. 

-  Geometric  scaling  often  has  no  sense,  e.g.  bubble  size  or 
surface  roughness  hardly  depend  on  system  dimensions. 

-  Of  the  proportion  problem,  arising  from  the  high  power 
density  levels,  being  typical  for  two  some  characteristic 
dimensions,  e.g.  bubble  size  and  two-phase  flow  and 
(boiling)  heat  transfer. 

The  approach  starts  by  summarising  the  results  of 
modelling  and  scaling  the  mechanically  or  capillary 
pumped  two-phase  loop,  shown  in  figure  1  (as  example). 


Figure  1.  Two-phase  loop  layout 


Conclusions  for  super-gravity  conditions  will  be  drawn 
from  this  summary  A  re-appraisal  of  the  dimensionless 
numbers  and  equations  has  to  identify  those  relevant  to 
model  and  scale  oscillating  heat  transfer  devices. 

DIMENSIONLESS  NUMBERS  &  SIMILARITY  ISSUES 

Similarity  considerations  (for  details  see  Refs.  1 ,  2)  led  to 
the  identification  of  1 8  dimensionless  groups  (so-called  7t- 
numbers),  considered  cruciai  for  the  thermal-gravitational 
scaling  of  two-phase  loops.  This  set  of  18  7t-numbers,  in 
literature  called  “the  complete  list  of  scaling  parameters” 
(Ref.  11),  is  shown  in  the  first  coiumn  of  Table  1.  The 
other  columns  show  the  different  sections  of  the  two- 
phase  heat  transport  system,  depicted  in  figure  1:  Pure 
liquid  lines  (with  and  without  heat  exchange),  capillary  or 
swirl  evaporators,  vapour  or  two-phase  lines,  and 
condensers.  For  each  section  it  is  indicated  by  •  whether  a 
typical  rc-number  is  relevant  in  that  section. 


Tablet. 

Relevance  of  n-numbers  for  thermal- 
gravitational  scaling  of  two-phase  loops 

Liquid  Parts 

Evaporators 
Swirl  & 
Capillary 

Non-liquid 

Lines 

Vapour/2-Phase 

Condensers 

Adiabatic 

Heating/Cooling 

m  =  D/L  =  geometry 

• 

• 

• 

« 

• 

m  =  Rei  =  (pvD/|j.)i  =  inertia/viscous 

• 

• 

• 

• 

• 

7t3  =  Fn  =  (v^/gD)i  =  inertia/gravity 

• 

• 

• 

/• 

• 

7t4  =  Eui  =  (Ap/pv^)i  =  pressure  head/inertia 

• 

♦ 

• 

• 

• 

^5  =  cos  V  =  orientation  with  respect  to  g 

• 

• 

• 

/• 

• 

JC6  =  S  =  slip  factor  =  Vv/vi 

• 

• 

• 

n?  =  density  ratio  =  pv/pi 

• 

• 

• 

Tts  =  viscosity  ratio  =  pv/pi 

• 

• 

• 

7C9  =  Wei  =  (pv*D/a)i  =  inertia/surface  tension 

• 

/• 

• 

n^o=  Pn  =  (pCpA,)i 

• 

• 

• 

jtii=  Nui  =  (hD/X.)i  =  convection/conduction 

• 

• 

• 

7Ci2=  Xv/k  =  thermal  conductivity  ratio 

• 

• 

jci3=  Cpv/Cpi  =  specific  heat  ratio 

• 

• 

ni4=  AH/hiv  =  Boil  ^enthalpy  nr.  =  X  =  quality 

• 

• 

• 

7115=  Moi  =  (pio^/pi''  g)  =  capillarity/buoyancy 

• 

/• 

• 

7t:ie=  Ma  =  v/(3p/3p)s’^^ 

• 

• 

• 

71)17=  {hAi)(m^g)''® 

• 

• 

7Ci8=  L^pi^  g  hivA,ipi(T-To) 

• 

• 
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There  is  perfect  similitude  between  model  and  prototype  if 
all  dimensionless  numbers  are  identical  in  prototype  and 
model.  Only  then  scaling  is  perfect.  It  is  evident  that  this  is 
not  possible  for  two-phase  flow  and  heat  transfer;  the 
phenomena  are  too  complex,  the  number  of  important 
parameters  or  rc-numbers  is  too  large.  Fortunately  also 
imperfect  (distorted)  scaling  can  give  useful  results  (Ref. 
12).  Therefore  a  careful  estimation  of  the  relative 
magnitudes  of  the  different  effects  is  required.  Effects, 
identified  to  be  of  minor  importance,  make  the  requirement 
for  identity  of  some  7i-numbers  superfluous  for  the  problem 
involved.  Examples  for  two-phase  systems  are;  the  Mach 
number  is  not  important  for  incompressible  liquid  flow,  the 
Froude  number  is  not  important  in  pure  vapour  flow. 

Further  it  can  be  remarked  that,  in  scaling  two-phase  heat 
transport  systems,  geometric  distortion  is  not  permitted  to 
study  boundary  layer  effects  and  boiling  heat  transfer,  as 
there  has  to  be  surface  roughness  identity  in  prototype 
and  model.  But  geometric  distortion  is  a  must  when  the 
length  scaling  leads  to  impractical  small  (capillary) 
conduits  in  the  model,  in  which  the  flow  phenomena 
basically  differ  from  flow  in  the  full  size  prototype. 

Sometimes  it  is  more  convenient  to  replace  quality  X  by  the 
volumetric  vapour  fraction  (void  fraction)  a,  according; 

(1 -a)/a  =  S(pv/pi)X/(1-X).  (1) 

It  is  clear  that  the  presented  set  of  7t-numbers  is  rather 
arbitrary,  e.g.  several  numbers  contain  only  liquid 
properties.  These  can  be  easily  transferred  into  vapour 
properties  containing  numbers  using  Tte  to  tcb.  Similarly 
can  be  used  to  interchange  a  characteristic  length  (duct 
length,  bend  curvature  radius)  and  some  characteristic 
diameter  (duct  diameter,  hydraulic  diameter,  but  also 
surface  roughness  or  bubble  diameter).  Sometimes  it  will 
even  be  convenient  to  simultaneously  consider  two 
geometric  tti -numbers.  One  concerns  the  overall  channel 
(channel  diameter  versus  length  or  bend  curvature  radius). 
The  second  pertains  to  other  parameters  as  the  ratio  of 
surface  roughness  and  bubble  diameter  to  investigate 
boiling  heat  transfer,  or  the  ratio  of  surface  roughness  and 
channel  diameter  to  study  friction  pressure  drop. 

The  best  scaling  approach  is  to  choose  (combinations  of) 
jt-numbers  such  that  they  optimally  suit  the  problem  under 
investigation,  e.g.; 

-  The  Morton  number  Mo, 

7ti5  =  Moi=RefFri/We®  =  pia®/p^g,  (2) 

which  is  useful  for  scaling  two-phase  flow  with  respect 
to  gravity,  as  it  contains  only  liquid  properties,  surface 
tension  and  gravity. 

-  The  Mach  number  Ma, 

jii6  =  Ma  =  v/(3p/ap)''",  (3) 

which  is  a  crucial  quantity,  if  the  compressibility  is 
important,  because  choking  depends  on  the  vapour 
quality  of  the  two-phase  mixture. 


-  The  enthalpy  or  boiling  number 

7ti4  =  Boil  =  AH(z)/  h|v=  Q/  tfihiv .  (4) 

Q  is  the  power  fed  to  the  boiling  liquid.  This  number 
appears  in  the  expression  for  the  dimensionless 
enthalpy  at  any  z  in  a  line  heated  from  outside  (q  is  the 
heat  flux); 

7ti4  =  AH(z)/h|v=  AHir/hiv  -H  TiDzq/rhhiv .  (5) 

For  sub-cooled  or  heated  liquid  this  is 

7Ii4  =  Q/rfi  Cpi  AT,  (6) 

AT  being  the  temperature  drop.  The  above  implies  that, 
if  the  dimensionless  entrance  enthalpies  are  equal  for 
different  fluids  flowing  in  a  similar  geometry,  equality  of 
the  boiling  number  ensures  equal  non-dimensional 
enthalpies  at  all  similar  axial  locations.  For 
thermodynamic  equilibrium  conditions  this  means  equal 
qualities  at  similar  locations,  and  similar  sub-cooling 
and  boiling  lengths. 

-  The  condensation  number,  in  which  h  is  the  local  heat 
transfer  coefficient, 

jci7  =  (h/A„(p,"/gp, (7) 

-  The  vertical  wall  condensation  number,  with  To  as  the 
local  sink,  T  as  the  local  saturation  temperature. 

71,8  =  L^Pi®  g  hiv/Mi  A,(T-  To) .  (8) 

A  first  step  in  a  practical  approach  to  scale  two-phase 
heat  transport  systems  is  identification  of  important 
phenomena,  to  obtain  7t-numbers  for  which  identity  in 
prototype  and  model  must  be  required  to  realise  perfect 
scaling  according  to  the  so-called  Buckingham  pi 
theorem  (crucial  in  similarity  considerations).  Distortion 
will  be  permitted  for  7t-numbers  pertaining  to  less 
important  phenomena.  Important  phenomena  and  the 
relevant  7t-numbers  will  be  different  in  different  parts  of  a 
system.  As  said  before,  the  relevance  of  the  7c-numbers 
in  the  various  loop  sections  is  indicated  by  •  in  Table  1  (tc- 
numbers  for  thermal-gravitational  scaling  of  two-phase 
loops. 

For  refrigerants,  like  ammonia  and  R114,  forced 
convection  heat  transfer  overrules  conduction  completely. 
Therefore  tiiq,  tch  and  7ti2,  are  not  critical  in  gravitational 
scaling.  Uis  can  be  neglected  also  as  the  system  maximum 
power  level  and  line  diameters  correspond  with  flow 
velocities  far  below  the  sonic  velocity  in  all  system  parts. 

Considering  713/7C5,  it  can  be  remarked  that  inertia  overrules 
buoyancy  not  only  in  pure  vapour  flow  or  in  a  low  gravity 
environment,  but  also  for  horizontal  liquid  sections  on  earth 
(v->7t/2).  This  implies  that  there  is  Tc-number  identity  for 
these  sections  in  the  low-g  prototype  and  in  the  terrestrial 
model,  for  a  horizontal  arrangement  of  these  sections. 
Also  it  is  remarked  that,  in  the  porous  (liquid)  part  of  a 
capillary-pumped  evaporator,  surface  tension  (a/Dp)  is 
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dominant  over  inertia  (719  ->  0):  hence  the  evaporator  exit 
quality  will  approach  1  (pure  vapour).  This  means  that 
gravity  is  less  important  for  the  vapour  part  of  the 
evaporator  and  the  vapour  line  to  the  condenser. 

A  first  conclusion,  that  can  be  drawn  now,  is:  Condensers 
and,  in  mechanically  pumped  systems,  also  two-phase 
lines,  are  crucial  in  scaling  with  respect  to  gravity.  They 
set  conditions  for  evaporators  and  single-phase  sections. 
A  second  conclusion  is:  In  adiabatic  two-phase  lines  of 
mechanically  pumped  systems  in  low-gravity,  only  shear 
forces  will  cause  separation  of  phases  in  a  high-quality 
mixture.  This  leads  to  annular  flow  (a  fast  moving  vapour  in 
the  core  and  a,  by  frictional  drag  induced,  slowly  moving 
liquid  annulus  at  the  inner  line  wall)  for  the  lower  flow  rates. 
For  increasing  power,  hence  flow  rate,  the  slip  factor  will 
increase  introducing  waves  on  the  liquid-vapour  interface 
and  entraining  of  liquid  droplets  in  the  vapour:  wavy- 
annular-mist  flow.  A  similar  flow  pattern  can  be  predicted 
for  vertical  downward  flow  on  earth,  as  it  easily  can  be 
derived  from  the  flow  pattern  map  for  downward  two-phase 
flow  shown  in  figure  2  (Ref.  13).  In  this  figure,  water 
properties  at  20  °C  must  be  used  to  determine  the  scale  of 
the  abscissa.  The  Froude  number  for  two-phase  flow  used 
in  this  figure  is  defined  as: 

Fr,p  =  (1 6  iti  VD®g)[X®/p^-i-  (1  -X)®/p®] .  (9) 


Figure  2.  Flow  pattern  map  (vertical  downward  flow) 

Comparing  low-g  and  vertical  downward  terrestrial  flow 
one  has  to  account  in  the  latter  for  the  reduction  of  the 
slip  factor  by  the  gravity  forces  assisting  the  liquid  layer 
flowing  down  (draining  effect).  Anyhow,  vertical  down  flow 
is  the  preferred  two-phase  line  orientation  in  the  terrestrial 
model,  because  of  its  axial-symmetric  flow  pattern. 


A  similar  conclusion  can  be  drawn  for  the  straight  tube 
condenser.  In  condensers  the  flow  will  change  from  wavy 
annular  mist  to  pure  liquid  flow,  passing  several  flow 
patterns,  depending  on  the  path  of  the  condensation. 

Consequences  of  scaling  are  elucidated  by  the  figures  3 
and  4,  showing  the  temperature  dependence  of  g.Moi  = 
Pl.a"/^i,‘'and  =  C.g*"  /(We/Fr)’^=  D.g’"  /(E6)’^=  C.g"" 
/2(E6)''^.  Eotvos  number  E6  and  Bond  number  Bo  are: 

E6  =  4Bo  =  gD^(Pl-Pv)/o.  (10) 


Temperature  (°C) 

Figure  3.  Pi.a^/pi^  versus  temperature  for  some  fluids 

The  conclusions  on  scaling  possibilities  were  drawn  from 
figures  like  the  figures  3  and  4  in  earlier  publications 
(Refs.  1  to  8).  They  pertain  to  scaling  at  the  same  gravity 
level  or  to  scaling  of  a  prototype  in  a  lower  than  1-g 
environment,  and  a  model  on  earth,  preferably  operating  in 
the  vertical  down-flow  (gravity  assist)  mode. 

They  can  be  summarised  by: 

-  Scaling  two-phase  heat  transport  systems  is  very 
complicated.  Only  distorted  scaling  offers  some 
possibilities,  when  not  the  entire  loop  but  only  loop 
sections  are  involved. 

-  Scaling  with  respect  to  gravity  is  hardly  discussed  in 
literature.  Some  possibilities  can  be  identified  for  very 
limited,  ti^ical  conditions  only. 
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-  Scaling  at  the  same  gravity  level  can  cover  only  a 
limited  range.  Scaling  of  high-pressure  system  (parts), 
e.g.  ammonia  at  110  °C,  by  low-pressure  system 
(parts),  e.g.  -50  °C  ammonia,  might  be  attractive 
either  for  safety  reasons  or  to  reduce  the  impact  of 
earth  gravity  in  vertical  two-phase  sections.  It  follows 
from  the  figures,  that  the  length  scale  ratio  between 
high-pressure  prototype  and  low-pressure  model 
(both  characterised  by  pi  a®/p'‘i  =  2.1  o’^  m/s^)  is  Lp/L^ 
=  [(a/pi)p/a/pi)m]’^^  =  0.4.  Anticipating  later  discussions, 
it  is  remarked  that  in  low  temperature  vertical  down 
flow  on  earth  the  impact  of  gravity  is  minimised, 
hence  it  is  suitable  to  simulate  conditions  below  1-g. 

-  Scaling  with  respect  to  gravity  is  restricted  to  say  two 
decades,  if  the  prototype  and  model  fluid  is  the  sahie. 

-  Scaling  a  full  size  low-gravity  (<  1 0'^  g)  mechanically- 

pumped  R1 14  loop,  like  the  ESA  mechanically 
pumped  Two-Phase  Heat  Transport  System,  can  be 
adequately  done  by  NLR’s  ammonia  test  rig,  since 
the  10'^-10'^g  R114  prototype  and  the  terrestrial 
ammonia  model  have  approximately  identical  Morton 
numbers.  This  fluid  to  fluid  scaling  leads  to  a  length 
scaling  Dp/Dm  =  (gm/Qp)  ^  =  4.5  to  6.5, 

in  agreement  with  the  ratio  of  actual  diameters:  21  mm 
for  the  R114  space  prototype  and  4.93  mm  for  the 
terrestrial  ammonia  model. 

-  NLR’s  mechanically  pumped  two-phase  ammonia 
test  rig  offers  also  some  opportunities  to  scale  the 
TPX  ammonia  loops  (Ref.  14). 

-  A  very  attractive  scaling  possibility  is  the  scaling  of  a 
two-phase  prototype  for  a  Mars  or  a  Moon  base,  by  a 
terrestrial  model  with  the  same  or  a  scaled  working 
fluid.  As  the  ratio  of  gravity  levels  between  prototype 
and  model  is  less  than  one  order  of  magnitude  (Mars 
0.4,  Moon  0.16),  the  sizes  of  the  model  have  to  be  only 
slightly  larger  than  the  geometric  sizes  of  the  prototype. 
In  addition,  adjustment  of  the  inclinations  (cos  v)  of 
non-horizontal  lines  in  the  terrestrial  model  may  help  to 
realise  almost  perfect  scaling. 

PRESSURE  DROP  &  HEAT  TRANSFER  EQUATIONS 


An  important  quantity  is  the  pressure  drop  in  the  different 
two-phase  sections  considered  crucial  for  two-phase  loop 
modelling  and  scaling.  The  next  considerations  therefore 
will  concentrate  on  pressure  drops  in  condensing  and 
adiabatic  flow  and  restrict  the  discussion  to  straight  tubes. 


The  total  local  pressure  gradient  for  annular  flow 
(dp(z)/dz)t  is  the  sum  of  friction,  momentum  and  gravity 
gradients.  Following  references  5  to  8,  and  15  (an 
elaborate  article  on  the  subject),  the  contribution  of 
friction  is  (deleting  the  z-dependence  to  shorten  the 
notation): 


(dp/dz),  =  -(32mVpvD“)(0.045/Rev''-^)[X'‘’+5.7(pi/iJv)' 


0.2\rx/1.8 . 


x0.0523^ 


*(1 -X)°''^X’ ®®(Pv/pi)°-^®^+8. 1  (pi/pv)°'^°®* 
H=(1-X)°®V®®(Pv/p|)°-®“].  (11) 


(12) 


Figure  4.  (a/p,  =  D.g''^.(We/Fr  versus  temperature, 

for  various  fluids 

Fluid  properties  are  assumed  independent  of  z,  since  they 
depend  only  on  the  mixture  temperature,  which  usually  is 
almost  constant  in  adiabatic  and  condensing  sections. 

The  momentum  constituent  can  be  written  as 

(dp/dz)m=  -(16rTf/7t^D''){[2X(1-a)/pva®  -  P(1-X)/pia+(1-P)* 

*(1  -X)/p|(1  -a)+(1  -X)/pi(1  -a)j(dX/dz)+ 

-  pC(1  -a)/pva®+(1  -X)^/pi(1  -a)^(da/dz)}.  (1 3) 

a  is  the  z-dependent  local  void  fraction  a(z).  p  =  2  for 
laminar  liquid  flow,  1 .25  for  turbulent  flow. 

The  gravity  constituent  is 

(dp/dz)g  =  ( 1  -a)(prPv)g  cosv.  (14) 

g  0  for  microgravity  conditions  and  g  cosv  equals  9.8 
m/s^  for  vertical  down  flow  on  Earth,  3.74  m/s^  for  vertical 
down  flow  on  Mars  and  1.62  m/s^  on  the  Moon,  a  is 
eliminated  in  (13)  and  (14)  by  inserting  (1). 

Slip  factor  S  is  to  be  specified.  The  principle  of  minimum 
entropy  production  (Ref.  16)  yields 

S=[(1+1.5Z)(p,/pv)]^''.  (15) 

This  is  for  annular  flow,  in  which  the  constant  Z  (according 
to  experiments)  is  above  1  and  below  2. 

S  =  {(Pi/Pv)[1  +Z'(pv/pi)(1  -X)/X]/[1  +Z'(1  -X)/X]}''^  (1 6) 

for  real  annular-mist  flow,  annular  flow  with  a  mass  fraction 
Z'  of  liquid  droplets  entrained  in  the  vapour.  Z'  is  between  0 
(zero  entrainment)  and  1  (full  entrainment).  In  the  limiting 
cases  Z  ->  0  and  Z'  0,  (15)  and  (16)  reduce  to 

S  =  (p,W. 


X  is  local  quality  X(z).  Re^  is  the  Reynolds  number 
Rev  =  4m/7cDpv. 


(17) 
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This  represents  ideal  annuiar  flow.  It  will  be  used  here  for 
simplicity  reasons  and  since  it  allows  comparison  with  the 
results  of  calculations  found  in  literature.  The  influence  of 
Z^O  and  Z'  9^  0  is  interesting  for  future  investigations. 

Inserting  (17)  into  (1)  and  (1 1, 13, 14),  yields 

(dp/dz),„  =  -(32itf/TT^PvD®).(D/2).(dX/dz)  [2(1-X)(Pv/p,)^®  + 
+2(2X-3+1/X)(Pv/p,r  +  (2X-1-|3X)(Pv/p,)^'"  + 

+  (2P-  pX-p/X)(Pv/pr®  +  2(1-X-P+pX)(Pv/p,)].  (18) 

(dp/dz)g  =  (32m2/TT^PvD®){1-I1+(Pv/pi)^"  (1-X)/X]''}* 

*[TT^D®g  cosv(Pi-Pv)Pv/32itP].  (19) 

To  solve  (11,  18,  19)  an  extra  relation  is  necessary, 
defining  the  z-dependence  of  X.  A  relation  often  used, 

dX/dZ  =  -Xen.rance/Lc.  (20) 

(Lc  =  condensation  length).  But  this  means  uniform  heat 
removal  (hence  linear  quality  decrease  along  the  duct), 
which  is  unrealistic.  It  is  better  to  use 

m  hiv(dX/dz)  =  -  h7tD[T(z)-Ts] ,  (21 ) 

relating  local  vapour  quality  and  heat  transfer,  h  is  the 
local  heat  transfer  coefficient  h(z),  for  which  we  can  write 

h  =  0.01 8(Aip/%i)Pri°®®l-(dp/dz),l’'®  (22) 

assuming  that  the  major  thermal  resistance  is  in  a  laminar 
sub-layer  of  the  turbulent  condensate  film  (Ref.  15). 

As  already  mentioned  the  two-phase  flow  path  is  almost 
isothermal,  which  implies  constant  temperature  drop  T(z)  - 
Ts  (for  a  constant  sink  temperature  Ts),  constant  fluid 
properties  and  a  constant  Prandtl  number,  defined  by 

Pr,  =  Cpi  P)  /A|  .  (23) 

The  total  condensation  pressure  drop  is 

U 

Apt=  J  (dp/dz)tdz.  (24) 

0 

The  equations  (1,  11,  18,  19,  21,  22)  can  be  combined. 
This  yields  an  implicit  non-linear  differential  equation  in  the 
variable  X(z),  which  can  be  rewritten  into  a  solvable 
standard  form  for  differential/algebraic  equations 

F(dX/dz,  X)  =  0.  (25) 

Figure  5  compares  the  pressure  gradient  constituents  for 
ammonia  vertical  down-flow  at  two  temperatures.  The 
curves  prove  that  the  impact  of  gravity  decreases  with 
decreasing  temperature.  This  confirms  the  earlier 
statement  that  low-gravity  behaviour  can  be  simulated  the 
best  by  terrestrial  down-flow  tests  at  low  temperature. 

Modelling  and  calculations  were  extended  from  adiabatic  to 
condensing  flow  in  a  straight  duct  (Refs.  17,  18),  in  order  to 
investigate  the  impact  of  gravity  level  on  the  duct  length 


required  to  achieve  complete  condensation.  This  impact, 
reported  to  lead  to  duct  lengths  being  more  than  one  order 
of  magnitude  larger  for  zero  gravity,  as  compared  to 
horizontal  orientation  on  earth  (Ref.  1 9),  was  assessed  for 
various  mass  flow  rates,  duct  diameters  and  thermal 
(loading)  conditions,  for  ammonia  and  R114.  A  summary 
of  results  of  calculations  for  ammonia  is  presented  next. 
To  compare  the  results  of  calculations  with  data  from 
literature,  the  condenser  in  reference  19,  was  chosen  as 
the  baseline.  Main  characteristics  are:  power  1  kW,  line 
diameter  16.1  mm,  ammonia  temperature  300  K  and 
temperature  drop  to  sink  10  K.  The  gravity  levels 
considered  are:  zero  gravity  g=0,  Earth  gravity  (1  -g)  g=9.8 
m/s^.  Mars  gravity  g=3.74  m/s^.  Moon  gravity  g=1 .62  m/s^ 
and  2-g  super-gravity  19.6  m/s^.  Illustrative  results  of 
calculations  are  discussed  next. 


Local  Vapjour  Quality  X 

Figure  5.  Pressure  gradient  constituents  as  function 
of  vapour  quality,  for  ammonia  at  +25‘C  and  -25*C 

Figure  6  depicts  the  vapour  quality  X  along  the 
condensation  path  (as  a  function  of  non-dimensional 
length  z/D)  for  all  gravity  levels  mentioned,  including  the 
curves  for  zero-g  and  horizontal  condensation  on  earth, 
found  in  literature  (Ref.  19).  From  this  figure  it  can  be 
concluded  that:  the  length  required  for  full  condensation 
strongly  increases  with  decreasing  gravity.  Zero-gravity 
condensation  length  is  roughly  10  times  the  terrestrial 
condensation  length.  The  data  of  reference  19  clearly  can 
be  considered  as  extremes. 

To  assess  the  impact  of  the  saturation  temperature  on 
condensation,  similar  cun/es  were  calculated  for  two  other 
temperatures,  243  K  and  333  K,  and  the  above  parameter 
values.  They  indicate  that  the  full  condensation  length 
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increases  with  tennperature  for  zero-g  conditions,  but 
decreases  with  temperature  for  other  gravity  ievels.  This 
impiies  that  differences  between  earth  gravity  and  iow-g 
decrease  with  decreasing  temperature,  it  confirms  the 
remark  that  gravity  impact  is  smaller  at  lower  temperatures. 


Nondimensional  length  z/D 


Figure  6.  Vapour  quality  along  the  reference  duct 

Calculations  of  the  vapour  quality  distribution  along  the  16.1 
mm  reference  duct  for  condensing  ammonia  (at  300  K) 
under  Earth  gravity  and  0-g  conditions,  for  power  levels 
ranging  from  0.5  kW  up  to  25  kW,  yielded  (Ref.  17)  that: 

-  A  factor  50  in  power,  25  kW  down  to  500  W, 
corresponds  in  a  zero  gravity  environment  to  a  relatively 
minor  reduction  in  full  condensation  length,  i.e.  from  600 
D  to  400  D  (9.5  to  6.5  m). 

-  Under  earth  gravity  conditions,  power  and  full 
condensation  length  are  strongly  interrelated:  from  U  = 
554  D  at  25  kW  to  only  1 9  D  at  500  W. 

-  The  gravity  dependence  of  the  full  condensation  length 
decreases  with  increasing  power,  until  the  differences 
vanish  at  roughly  1  MW  condenser  choking  conditions. 
The  latter  value  is  an  upper  limit,  calculated  (Ref.  16)  for 
ideal  annular  flow.  Choking  may  occur  at  considerably 
lower  power  values  in  the  case  of  actual  annular-wavy- 
mist  flow,  but  the  value  exceeds  anyhow  the  choking 
limit  for  homogeneous  flow,  roughly  170  kW. 

Calculation  of  the  vapour  quality  along  the  duct  for  three 
gravity  levels  (0,  Earth  and  2-g)  and  three  duct  diameters 
(8.05,  16.1  and  24.15  mm)  at  300  K,  yielded  the  ratio  of  the 
absolute  duct  lengths  Lc(m)  needed  for  full  condensation 
under  zero-g  and  one-g  respectively  (Ref.  17).  It  has  been 
concluded  that  the  ratio  between  full  condensation  lengths 
in  zero-g  and  on  Earth  ranges  from  roughly  1 .5  for  the  8.05 
mm  duct,  via  1 1  for  the  16.1  mm  duct,  up  to  more  than  30 
for  the  24.15  mm  duct.  In  other  words,  small  line  diameter 
systems  are  less  sensitive  for  differences  in  gravity  levels 
as  compared  to  larger  diameter  systems.  This  has  been 
confirmed  by  TPX  I  flight  data  (Ref.  20). 


As  the  model  developed  is  mainly  valid  for  annular  flow,  it 
will  be  worthwhile  to  investigate  the  impact  of  other  flow 
patterns  inside  the  condenser  duct  (mist  flow  at  high 
quality,  slug  and  bubbly  flow  at  low  quality  and  wavy- 
annular-mist  in  between).  In  other  words,  it  is  to 
investigate  whether  the  pure  annular  flow  assumption, 
leads  towards  slightly  or  substantially  overestimated  full 
condensation  lengths.  Complications  are  that  the 
boundaries  flow  patterns  are  not  accurately  known  and 
flow  pattern  transitions  strongly  depend  on  temperature 
and  line  diameter.  Accurate  knowledge  of  the  gravity  level 
dependent  two-phase  flow  regimes  is  crucial  for  modelling 
and  designing  two-phase  heat  transport  systems  for 
space,  as  flow  patterns  directly  affect  thermal  hydraulic 
characteristics  of  two-phase  flow  and  heat  transfer. 
Therefore  flow  pattern  (regime)  maps  are  to  be  created, 
preferably  in  the  non-dimensional  format  (Fig.  2)  or  in  the 
3-dimensional  format  of  the  figures  7  and  8,  presenting 
data  measured  with  a  R12,  10.5  mm  line  diameter  two- 
phase  loop  during  many  K135  aircraft  flights  (Ref.  21). 
Such  flow  pattern  maps  can  then  be  used  to  determine,  in 
an  iterative  way,  via  the  flow  pattern  dependent  constitutive 
equations  for  two-phase  flow  and  heat  transfer,  the  actual 
trajectories  of  condensing  or  evaporating  (boiling)  flow. 
The  latter  will  finally  lead  to  an  accurate  determination  of 
the  pressure  drops  in  the  various  sections  and  of  the  heat 
transfer  in  the  evaporator  or  condenser  sections  of  a  two- 
phase  heat  transport  system. 


Figure  7.  Annular  flow:  Gravity  dependent  3-D  map 


Figure  8.  Slug/plug  flow:  Gravity  dependent  3-D  map 


-10- 

NLR-TP-2000-213 


Figure  9  qualitatively  illustrates  the  aspects  of  pressure 
gradient  characteristics  of  various  flow  regimes  in  terms 
of  an  apparent  interfacial  friction  factor. 


Figure  9.  Qualitative  aspects  of  pressure  gradient 
(apparent  friction  factor)  for  various  flow  patterns 

In  summary  it  can  be  said  that  the  information  presented 
confirms  the  results  of  other  models.  When  designing 
condensers  for  space  applications,  one  should  carefully 
use  and  interpret  data  obtained  from  terrestrial  condenser 
tests,  even  when  the  latter  pertain  to  vertical  downward 
flow  situations  (characterised  by  the  same  flow  pattern). 

The  equations  given  are  useful  for  a  better  understanding 
of  the  problems  that  can  be  expected:  e.g.  the  needed 
condenser  lengths  in  space  applications.  Equations  and 
results  of  the  calculations  suggest  that  hybrid  scaling, 
which  combines  geometric  and  fluid-to-fluid  scaling,  can 
usefully  support  the  design  of  space-oriented  two-phase 
heat  transport  systems  and  their  components.  With  respect 
to  the  local  heat  transfer  equation  used  in  equation  (22),  it 
is  remarked  that  it  has  a  wrong  lower  limit  h^O  for 
(dp/dz)t^0,  which  disappears  by  incorporating  conduction 
via  the  liquid  layer.  Preliminary  calculations  indicate  that 
incorporation  of  pure  conduction  will  lead  to  somewhat 
shorter  full  condensation  lengths,  both  for  zero  and  for  non¬ 
zero  gravity  environments.  This  implies  only  quantitative 
changes:  Consequently  the  conclusions  presented  before 
remain  valid. 

Finally  it  is  stressed,  that  the  design  of  ammonia  two-phase 
heat  transport  system  prototypes  for  Mars  and  Moon  base 
applications,  looks  very  promising.  This  is  because  flow 
pattern  maps  for  planetary  gravity  levels  can  be  easily 
obtained  from  measurements  in  models  on  earth,  with  the 
same  or  with  another  working  fluid.  These  models  will 
have  (almost)  identical  geometry  and  (if  necessary)  equal 
lines.  They  also  might  have  a  dedicated  inclination  with 
respect  to  the  Earth  gravity  vector  (i.e.  cos  v  =  0.4  in  the 
case  of  Mars,  0.16  in  the  case  of  the  Moon).  These  maps 
can,  in  addition,  include  data  obtained  from  experiments 
with  ammonia  test  loops,  carried  out  during  low-g  aircraft 
trajectories. 


1-G  SCALING  OF  PROTOTYPES  FOR  SUPER-G 

The  preceding  considerations,  pertaining  to  gravity-assist 
conditions  (vertical  down-flow  in  some  gravity  field),  remain 
valid  and  useable  for  two-phase  loops,  including  pulsating 
two-phase  loops,  in  super-gravity-assist  conditions. 

A  quantitative  example  is  the  simple  extrapolation  of  the 
gravity  dependent  condensation  curves  (Fig.  6)  to  12-g, 
which  yields  a  1 0-g  gravity-assist  fuli  condensation  length 
of  the  order  of  10  D. 

Though  many  things  will  be  different  in  the  anti-gravity 
mode,  the  general  thermal-gravitational  scaling  ruies  are 
valid  for  both  gravity  assist  and  anti-gravity  conditions. 
Two  interesting  scaling  possibilities,  directly  derivable 
from  the  figures  3  and  4,  are  discussed  hereafter. 

First,  ammonia  is  the  most  promising  working  fluid  for  two- 
phase  systems  operating  between  270  and  350  K,  e.g.  for 
cooling  electronics  in  a  super-gravity  environment  around 
say  10-g.  Water,  having  better  performance  above  350  K, 
is  also  a  candidate  fluid.  But  it  can  freeze  when  the  system 
is  not  operating  (which  is  undesirable  or  unacceptable). 
Also  its  performance,  being  under  315  K  below  the 
acetone  performance,  becomes  even  unacceptably  low 
under  300  K. 

Figure  3  yields  for  ammonia  in  the  aforementioned 
operation  range  (~  10-g  at  270  -  350  K)  a  value  of  g.Moj 
between  1  - 1.4*10^®.  This  implies  that  the  working  fluid  of 
the  terrestrial  model  must  have  a  g.Moi  of  1  -  1.4*10^^. 
This  requirement  is  fulfilled  by  acetone,  within  roughly  the 
same  temperature  range.  Figure  4  straightfon/vardly  gives 
that  the  corresponding  geometric  scaling  ratio  between 
model  and  prototype  (Lm/Lp  =  Dm/Dp)  is  between  3  and  3.5. 
The  above  features  make  acetone  well  suited  for  1-g  scale 
models  to  simulate  10-g  ammonia  prototypes 

Secondly,  if  ammonia  is  not  an  acceptable  working  fluid 
(e.g.  because  of  safety  requirements)  and  water  is  not  a 
candidate  (because  of  the  earlier  mentioned  reasons), 
acetone  might  be  the  alternative  for  a  10-g,  270  -  350  K 
prototype.  Figure  3  yields  that  the  g.Moi  of  1  -  1.4*10^^  for 
acetone  can  be  adequately  scaled  by  R114  at  1-g  within 
the  same  temperature  range  or  a  slightly  lower  one.  The 
corresponding  geometry  ratio  between  model  and 
prototi^e  (Lm/Lp  =  Dm/Dp)  is  about  1 .5,  according  to  figure 
4.  Consequently,  R114  can  be  considered  useful  for  1-g 
scale  models  to  simulate  10-g  acetone  prototypes. 

ANTI-GRAVITY  ISSUES 

The  impact  of  anti  or  against  gravity  or  super-gravity  can 
be  assessed  by  recalling  figure  6,  showing  the  absolute 
values  of  the  three  pressure  drop  constituents.  The  results 
of  pressure  drop  calculations,  presented  in  the  preceding 
sections,  pertained  to  gravity-assist  condensation  This 
implies  that  the  calculations  were  done  for  cosv  =  +1, 
meaning  (as  it  is  gravity-assist)  that  the  contribution  of 
gravity  had  the  opposite  sign  of  the  constituents  of  friction 
and  momentum. 
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The  results,  of  the  only  calculations  for  a  case  where  the 
gravitational  acceleration  was  against  the  flow  direction  of 
the  two-phase  mixture,  are  depicted  in  figure  10. 


Quaiity  X 

Figure  10.  Measured  versus  calculated  adiabatic 
pressure  drops  in  a  R114  duct 

The  figure  depicts  the  flight  data  (Ref.  22)  and  three 
curves:  one  with  the  acceleration  in  the  flow  direction,  one 
against  the  flow  direction,  and  one  perpendicular  to  the 
flow.  The  latter  were  considered  to  be  extremes,  as  only 
the  gravity  level  (10'^g)  was  published,  not  the  direction. 
Though  10  ®g  is  only  a  small  acceleration,  the  curves 
clearly  illustrate  an  impact  of  reversai  of  the  direction  of 
gravity.  It  is  obvious  that  this  impact  will  be  orders  of 
magnitude  larger  at  super-gravity  leveis. 

The  above  suggests  that  useful  anti-gravity  performance 
calculations  can  be  carried  out.  But  one  shall  keep  in 
mind  that  the  equations  used  pertain  to  annular  flow  only. 

When  considering  anti-(super)-gravity  modes,  four 
important  issues  are  to  be  noted: 

-  First,  any  pump  of  a  two-phase  system  has  to  deiiver 
a  pumping  pressure  that  is  at  ieast  the  sum  of  ail 
pressure  losses  in  the  systems,  in  an  anti-1  Og  mode, 
this  is  mainly  the  gravitational  pressure  drop.  Figure  5 
illustrates  that  the  contributions  of  friction  and 
momentum  become  negligible,  as  compared  to  the 
12-g  contribution,  which  can  by  obtained  by  simply 
shifting  the  1-g  curve  one  decade  upwards.  The  10-g 
curve  clearly  overwhelms  the  curves  of  the  other 
constituents,  for  almost  all  vapour  qualities  except  for 
quaiities  over  say  0.9,  where  the  flow  pattern  probably 
will  be  homogeneous,  instead  annular.  In  case  of  a 
mechanical  pump,  the  huge  (super-gravity)  pressure 
losses  can  be  overcome  by  the  use  of  a  correctly 
designed  displacement  pump.  Existing  capillary 
pumped  systems,  however,  will  fail  completely,  since 
existing  capillary  structures  can  operate  a  two-phase 
loop  against  gravity,  only  if  the  maximum  pumping 


height  to  be  delivered  is  less  than  say  6  metres.  The 
6  m  gravitational  pressure  drop  is  at  least  a  factor  20 
smaller  than  the  10-g  super-gravity  pressure  drops. 
Current  high-performance  capillary  pumps  use 
reasonably  permeable  wicks  with  a  pore  size  of  0.5 
pm  to  realise  the  6  metres  pumping  height.  For  a  20 
times  larger  pumping  height,  the  pore  size  shall  be 
about  25  nm.  But  a  wick  of  such  small  pores  will  not 
have  the  required  permeability  of  a  capable  capillary 
pump.  Consequences  for  vapour  pressure  driven 
two-phase  systems,  iike  oscillating  (pulsating)  heat 
transport  devices,  wili  be  discussed  in  a  later  section. 

-  Second,  the  pressure  and  temperature  drops  are 
only  minor  in  the  gravity-assist  condensation  curves 
presented.  It  means  that  the  condensation  is  almost 
isothermal,  hence  the  fluid  properties  can  be 
assumed  constant.  This  will  not  be  the  case  for  anti- 
super-gravity  conditions,  meaning  that  for  the  latter 
conditions  calculations  will  become  very  complicated. 

-  Third,  for  super-g  conditions  there  is  no  information 
at  all  on  flow  pattern  maps  and  the  boundaries 
between  the  different  flow  regimes.  Therefore  and  as 
it  is  expected  that  the  various  items  will  substantially 
differ  from  the  (hardly  available)  existing  1-g  ones, 
the  creation  of  consistent  flow  pattern  maps  for  the 
super-gravity  environment  has  been  started  at  NLR. 

-  Fourth,  gravity  influences  the  pressure  drop  and  the 
corresponding  temperature  drop  across  the  heat 
transport  system.  This  is  confirmed  by  experimental 
data  (Ref.  9):  Compared  to  0-g  pressure  and 
temperature  drops,  the  drops  decrease  with 
increasing  gravity-assist,  and  increase  with 
increasing  anti-gravity. 

SINGLE-PHASE  OSCILLATING  HEAT  TRANSFER 

The  assessment  of  super-gravity  performance  issues  of 
osciilating  (pulsating)  heat  transfer  devices  starts  with  the 
description  of  a  very  interesting  device  (Fig.  11),  foilowing 
a  short  earlier  presented  summary  (Ref.  23). 

For  detailed  discussions  on  this  synchronised  forced 
oscillatory  flow  heat  transfer  device,  it  is  referred  to  the 
publications  of  the  originators  (Refs.  24-26)  and  to  other 
publications  on  this  device  and  on  the  related  phase- 
shifted  forced  oscillatory  flow  heat  transfer  devices  (Ref. 
27). 

The  device  is  not  a  two-phase  device  but  a  single-phase 
one.  The  set-up  consists  of  two  reservoirs  at  different 
temperatures,  connected  by  a  0.2  m  long,  1 2.7  mm  inner 
diameter  acryiic  tube,  containing  31  glass  capillaries  with 
an  inner  diameter  d  =  1  mm.  The  open  cross-sectional 
area  of  the  capillary  structure,  including  the  triangular 
sections  between  the  capillaries,  Ai  was  determined  to  be 
67  mm^,  being  53%  of  the  total  inner  cross-sectional  area 
of  the  tube  (A  =  127  mm^*  The  reservoirs  are  equipped 
with  flexible  membranes.  A  variable  frequency  shaker  is 
used  to  oscillate  the  incompressible  working  liquid  inside 
the  capillary  structure.  The  frequency  f  is  variable  from  2 
and  8  Hz.  The  tidal  displacement  Az  is  variable  from  20 
to  125  mm. 
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The  operation  is  as  follows.  Starting  with  the  capiilary 
structure  fiiled  with  hot  iiquid,  this  iiquid  is  replaced,  in  the 
first  haif  of  the  osciilation  period,  by  liquid  from  the  coid 
reservoir,  but  with  the  exception  of  the  thin  (Stokes) 
boundary  iayer.  Heat  is  then  exchanged  very  effectiveiy  in 
radiai  direction  between  the  hot  Stokes  layer  and  the  cold 
core.  Heat  accepted  by  the  core  is  removed  to  the  cold 
resen/oir  in  the  second  half  of  the  oscillation  period.  The 
heat  flow,  via  the  liquid,  between  the  reservoirs  equals 


Q  =  Aeff  (A|/L)  AT  =  Pi  Cp,  KeH  (A/L)  AT.  (26) 

AT  is  the  temperature  difference  between  the  hot  and 
cold  reservoir.  Kew  is  the  effective  thermal  diffusivity ,  Aeft  is 
the  effective  thermal  conductivity. 


Flexible 

membrane 
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Fiexible 

membrane 


Figure  11.  Synchronised  forced  oscillatory  flow  heat 
transfer  device 

Figure  12  shows  the  experimentally  determined  effective 
thermal  diffusivity  as  a  function  of  tidal  displacement  and 
oscillation  frequency,  for  a  device  with  glass  capillaries 
and  water  as  working  fluid.  The  solid  lines  are  analytical 
predictions  from  the  laminar  theory.  The  figure  indicates 
that  the  effective  thermal  conductivity  via  the  liquid  is 

Aeff  =  B  Pi  Cp,  {{Azf  /  (d/2)}  {2k  f  p, /p,)''".  (27) 


Figure  12.  Effective  thermal  diffusivity 

An  important  conclusion  for  future  design  activities  can 
be  drawn  by  considering  the  equations  (28)  to  (31):  The 
highest  values  of  proportionality  factor  B  occur  for  small 
Prandtl  number  fluids  and  walls  which  are  good  thermal 
conductors  (Ref.  24). 

An  alternative  representation  is  obtained  by  defining  an 
enhancement  (proportionality)  factor  E  for  undeveloped 


oscillating  flow  in  synchronised  systems 
(Aefl  /  A|)  -1  =  (Keff  /  Ki)  -1  =  (Pr,  Az  /2d)".E  .  (31 ) 

E  depends  on  the  dimensionless  Womersley  number 
W0|  =  (d/2)(27cfpi/pi)''".  (32) 

It  can  be  approximated  by 

E  =  W0|''/24  for  Wo, « 1  ,  (33) 

E  =  W0|  (Wo,  -  2'^'^ )  for  Wo,  «  1  .  (34) 


Figure  13  shows  the  enhancement  factor  E  as  a  function 
of  Womersley  number,  for  decades  of  the  Prandtl  number 
0.1, 1  and  10  (Refs.  25, 28). 


Proportionality  factor  B,  the  tangent  of  the  straight  lines  in 
figure  1 0,  can  be  written  as  (Refs.  24,  26) 

B  =  2'®®  Pr,’’  B’  +  B”  -  B’B”  (1  +  Pr,^'®)  * 

*{Pri'^'^-2(1 +Pri)})i-.  (28) 

B’  =  Pri/(Pr|-1),  (29) 

B”  =  {(1  -  B’)  (Pr,  K1  /  Kwal,)"'^-  B’  A,  /  Awall)}  * 

*{A|/Awa„  +  (iq/K„all)'V-  (30) 


Inserting  properties  of  water  at  20  °C,  d  =  1  mm,  Az  = 
125  mm  and  f  =  8  Hz  (hence  Pr,  =  6.9  and  Wo,  =  3.55), 
yields  (according  to  figure  13)  for  the  above  synchronised 
device  an  enhancement  factor  0.14.  The  corresponding 
Aeff  (about  1 .2  *  1 0'*  W/m.K)  means  a  power  density  slightly 
above  3*10®  W/m^.  This  is  confirmed  by  the  experimental 
data:  2.9*10®  W/m^  for  a  gradient  of  280  K/m  (Ref.  24). 
The  total  effective  thermal  conductivity  is  roughly  5000 
W/m.K,  if  it  is  based  on  the  total  cross-sectional  tube 
area,  including  the  insulating  glass  walls. 


-13- 

NLR-TP-2000-213 


Figure  13.  Enhancement  factor  E  as  a  function  of 
Womersley  number  Woj  and  Prandtl  number  Pn 

A  similar  expression  has  been  derived  for  undeveloped 
oscillating  flow  in  phase-shifted  systems  (Refs.  27, 29): 

Aeff  =  A|  {1  +  0.707  (1  +  PrrV  (1  +  * 

*  {(Az)®/(d/2)}{(27C  f  Pi  CpA)'''}  .  (35) 

Inserting  the  properties  of  water  at  20  °C,  d  =  1  mm,  Az  = 
125  mm  and  f  =  8  Hz  yields  Aeff  =  2.5*10'*  W/m.K,  being 
2.5  times  the  value  of  the  corresponding  synchronised 
forced  oscillatory  flow  heat  transfer  device. 

Disadvantages  of  the  concepts  are  the  poor  current  state 
of  the  art  and  the  power  consumption  of  the  shaker.  The 
latter,  estimated  to  be  more  than  5  W  in  the  described 
cases,  strongly  increases  with  oscillation  frequency  and 
tidal  displacement.  A  disadvantage  for  some  micro-gravity 
applications  can  be  the  noise  introduced  by  the  shaker. 

The  major  advantage  of  the  concepts  is  its  variable 
conductance,  which  is  adjustable  via  the  frequency  and 
tidal  displacement  from  almost  zero  (liquid  in  rest  in  a 
poorly  conductance  structure)  to  values  comparable  to  or 
even  better  than  heat  pipes.  This  makes  such  devices 
very  useful  for  instance  to  drain  huge  amounts  of  thermal 
power  from  a  hot  vessel  for  nuclear  reactor  cooling,  in 
case  of  an  emergency. 

OSCILLATING  TWO-PHASE  HEAT  TRANSFER 

Two  oscillating  two-phase  heat  transfer  devices  can  be 
distinguished:  Pulsating  two-phase  heat  transfer  loops  and 
devices  like  the  pulsating  or  meandering  heat  pipe  and  the 
flat  swinging  heat  pipe.  These  systems  have  in  common 
that  the  operation  is  driven  only  by  the  vapour  pressure 
differences  induced  by  the  heat  to  be  transported  itself.  In 
other  words,  they  do  not  need  an  additional  power  source. 


Figure  14  depicts  the  temperature  dependent  saturation 
pressure  of  some  candidate  working  fluids.  It  illustrates 
that  when  designing  a  device,  one  has  to  select  a  working 
fluid  with  a  high  saturation  pressure  gradient  (dp/dT)  within 
the  operating  (temperature)  range,  as  higher  system 
pumping  pressures  correspond  with  higher  dp/dT-  values. 


Figure  14.  Vapour  pressures  of  candidate  working 
fluids,  as  a  function  of  the  temperature 


PULSATING  TWO-PHASE  LOOPS 

Pulsating  two-phase  heat  transfer  loops  already  were 
proposed  many  years  ago  (Refs.  29  to  31).  These  systems 
are  just  two-phase  loops  driven  by  vapour  pressure 
differences,  instead  of  mechanical  or  capillary  pumping 
action.  The  vapour  pressure  pumping  action  is  realised  by 
incorporation  of  two  one-way  valves:  one  at  the  entrance, 
the  other  at  the  exit,  of  a  separate  unit,  an  evaporator  or  a 
control  reservoir.  Power  fed  to  the  system,  increases  the 
internal  vapour  pressure  in  the  section  between  the  two 
valves  till  the  exit  valve  will  open  and  the  loop  starts  to  run, 
also  opening  the  second  valve  to  let  sub-cooled  liquid  flow 
into  the  pumping  section.  After  some  pressure  decay  the 
valves  will  close  re-starting  the  process.  A  careful  design 
will  certainly  lead  to  a  properly  performing  heat  transfer 
device.  The  advantages  of  such  loops  are  the  already 
mentioned  driving  mechanism  (the  heat  to  be  transferred 
only,  without  an  additionai  power),  the  high  heat  transport 
capability,  their  self-priming  capability,  and  their  capabiiity 
to  work  against  gravity.  A  disadvantage  is  the  pulsating 
operation,  hence  pulsating  heat  and  mass  transfer, 
accompanied  by  temperature  variations  and  possibiy  also 
vibrations  (g-jitter),  which  will  make  the  system  not 
attractive  for  some  micro-gravity  appiications. 

The  majority  of  the  remarks  on  two-phase  thermal- 
gravitational  modelling  and  scaiing  made  up  to  now  are 
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directly  applicable  for  pulsating  two-phase  loops.  The  only 
new  issue  to  be  accounted  for  is  the  proper  design  of  the 
desired  driving  pressure,  which  is  directly  coupled  to  the 
corresponding  difference  of  the  saturation  temperatures  in 
evaporator  and  condenser. 

OTHER  OSCILLATING  TWO-PHASE  DEVICES 

Other  oscillating  devices  have  various  names,  e.g.  bubble- 
driven  heat  transfer  device  (Ref.  27),  pulsating  or 
meandering  heat  pipe  (Refs.  32,  33),  capillary  heat  pipe  or 
capillary  tunnel  heat  pipe  (Refs.  34,  35),  flat  swinging  heat 
pipe  (Ref.  36),  and  spirally  wound  or  serpentine-like  heat 
pipe  (Ref.  37). 


hot  evaporator  side 


vapour 
'  bubble 


liquid 

'slug 


cold  condenser  side 

Figure  14.  Schematic  of  an  osciliating  device  section 


Figure  14  schematically  shows  a  section  of  such  an 
oscillating  device,  called  looped  (Ref.  38)  or  closed-loop 
(Ref.  32),  if  the  two  legs  at  each  end  are  not  dead  ends  but 
interconnected,  thus  creating  a  closed  loop  configuration.  If 
the  latter  configuration  has  a  spring-like  geometry  like  the 
arrangement  discussed  in  references  38  and  39,  the 
operation  of  the  device  has  been  frequently  observed  to 
stabilise  after  a  certain  start-up  period,  producing  a  periodic 
pumping  of  almost  constant  frequency  into  one  direction 
(Ref.  39).  Such  behaviour  is  very  similar  to  the  behaviour  of 
the  just  described  pulsating  heat  transfer  loops.  The 
similarity  suggest  that  the  function  of  the  valves  in  the 
pulsating  heat  transfer  loops  is  being  delivered  by  the  stick- 
slip  conditions  of  the  slug-plug  distribution  of  the  working 
fluid  inside  the  closed-loop  spring-like  meandering 
structure,  whose  heating  and  cooling  sections  have  also  a 
certain  periodicity.  Anyhow,  the  slug-plug  distribution  is 
essential  for  these  oscillating  devices,  without  valves. 


Consequently,  the  thermal-gravitational  scaling  of  the  inner 
tube  diameter  can  be  derived  from  figure  4. 

Second,  the  slug-plug  condition  also  sets  the  lower  limit  of 
the  slug  (bubble)  size:  It  shall  be  at  least  0.6  times  the  tube 
diameter  (Ref.  38).  This  extra  requirement  has  an  impact 
on  the  liquid  filling  ratio  (1  -  a).  If  not  fulfilled,  the  bubbles 
will  be  too  small  to  maintain  the  slug  flow  pattern, 
characterised  by  high  heat  transfer  densities.  The  resulting 
bubbly  flow  leads  to  a  far  less  efficient  heat  transport. 

Other  useful  plots  of  dimensionless  numbers  are  shown 
in  the  figures  16  and  17  (Ref.  39).  Figure  16  shows  the 
dimensionless  velocity  v*  as  a  function  of  the  Morton 
number  Mo  and  the  Eotvos  number  E6. 


Eo  =  (pjj^gD2/a) 


Figure  16.  Dimensionless  velocity  v*  as  a  function  of 
Morton  number  Mo  and  Eotvos  number  Ed 

Figure  17  depicts  many  experimental  data,  but  in  the 
alternative  plotting  (Ref.  40)  v*  as  a  function  of  the  inverse 
viscosity  number  Mu,  for  different  values  of  Archimedes 
number  Ar. 


The  velocities  of  bubbles  and  slugs  in  a  tube  are  governed 
by  buoyancy,  liquid  inertia,  liquid  viscosity  and  surface 
tension  forces,  in  the  general  case  of  bubbly  or  slug-plug 
flow  in  a  gravity  field.  This  means  that  properly  chosen 
dimensionless  groups  can  be  very  helpful  to  discuss  the 
aspects  of  slug-plug  flow  in  oscillating  devices  (Ref.  40). 

First  the  condition  of  slug-plug  distribution  determines  the 
maximum  inner  tube  diameter  (Refs.  40, 41 ,  32) 

Dr„ax  =  1 .836  g [o/(Pi -  Pv)]"'"  S  1 .836  g-"'"  (o/p,)''"  (36) 

forpi»pv. 


The  dimensionless  numbers  Mu  and  Ar  are  given  by 
Mu  =  p,  [g  (p,  -  Pv)  p,]  =  Pi  (g  D®  p,®)  (37) 

[krf  =  Mo  =  (Pi  o^/pi'*  g)  /  (1  -  Pv/Pil  s  p,  0®^  g  .  (38) 

The  three  asymptotes  shown  are: 

-  V*  =  0.345,  for  Ed  >  100  and  Mu  <  10■^  the  inertia 
dominant  domain. 

-  V*.  Mu  =  10'^,  for  Ed  >  100  and  Mu  >  0.5,  the  viscosity 
dominant  domain. 

-  Mu^  Ar  =  0.16  and  Ed  <  3.37,  the  surface  tension 
dominant  domain. 
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The  last  domain  is  the  most  important  for  the  oscillating 
devices  considered  here,  as  E6  <  3.37  straightforwardly 
leads  to  the  maximum  tube  diameter  equation  (30)  and 
dominating  surface  tension  means  that  piug-slugs  do  not 
move,  if  there  is  no  (thermal)  power  input.  In  this  domain 
slug  flow  is  guaranteed  by  surface  tension,  if  siug  bubbles 
have  diameters  of  at  least  0.6  times  the  tube  diameter. 


Figure  17  Dimensionless  velocity  v*  as  a  function  of 
dimensionless  inverse  viscosity  number  Mu 

Some  quantitative  considerations  on  these  pulsating  two- 
phase  heat  transfer  devices  are  presented  next,  in  order  to 
get  a  certain  feeling  for  or  a  better  understanding  of  how  to 
design,  scale  and  test  such  devices.  Geometrical  data  and 
some  performance  figures  of  the  earlier  described  single 
phase  oscillating  heat  transfer  device  will  be  used. 

EXPERIMENTAL  DATA  FROM  LITERATURE 

Several  recently  published  resuits  of  experiments  confirm 
the  remarks  made  up  to  now.  They  wili  be  summarised 
first,  before  being  used  to  establish  a  simple  quantitative 
model  for  oscillating  (pulsating)  heat  transfer  devices. 

Super-gravity  experiments  were  done,  on  a  4.5  m  diameter 
centrifuge  table,  with  a  non-looped  device,  consisting  of  23 
turns  of  0.42  m  long,  1.1  mm  internal  diameter,  stainless 
steel  capillaries  (Ref.  9).  The  working  fluid  was  acetone, 
the  liquid  filling  ratio  (1-a)  was  0.6.  The  length  of  the 
sections  was  120  mm  for  the  heating  and  the  adiabatic 
section,  180  mm  for  the  cooling  section.  The  experiments 
were  done  in  both  gravity-assist  and  anti-gravity  conditions. 
Experimental  data  confirms  that  gravity  influences  the 
pressure  drop  and  the  corresponding  temperature  drop 
across  such  devices.  For  example,  while  continuously 
transporting  the  maximum  power  40  W,  the  finally 
reached  stable  heater  section  temperature  increased 
from  130  °C  in  the  6-g-assist  case,  via  160  °C  at  0-g  and 
185  °C  at  6-g  anti-gravity,  to  200  °C  at  12-g  anti-gravity. 


There  occurred  no  evaporator  dry-out  for  all  acceleration 
conditions  specified.  The  above  tests,  and  vibration  tests 
(frequency  0  to  16  kHz,  acceleration  0  to  15  g,  amplitude  0 
to  7  mm,  inclination  0  - 1 80°),  showed  that  these  oscillating 
devices  are  not  sensitive  for  acceleration  fields. 

Experiments  with  a  4-turn,  396  mm  long,  glass,  closed- 
loop  meandering  device  showed  good  performance  also 
(Ref.  27),  even  for  capillary  diameters  larger  than  the 
maximum  diameter  D^ax  (see  equation  36).  But  poor 
performance  was  observed  for  diameters  below  0.6  Dmax. 

Results  of  experiments  with  a  2.4  mm  diameter  tube  test 
device  showed  for  all  liquids  a  maximum  performance  at  a 
liquid  filling  ratio  (1-a)  around  0.35.  For  filling  ratios  up  to 
0.8,  the  performances  are  slightly  below  these  maximum 
performances,  above  0.8  they  strongly  decrease.  Typical 
values  of  the  performance  are  about  1.25  W/K  for  water 
and  0.8  W/K  for  ethanol  and  R141b,  for  a  temperature 
drop  of  50  K,  within  the  working  range  10  -  100  °C.  The 
total  wall  resistance  of  the  heating  and  cooling  sections 
were  calculated  to  be  1 .25  K/W  per  tube. 

A  typical  periodic  vapour-plug  propagation  phenomenon 
was  reported  to  occur  in  a  10-tum  water  system  (Ref.  32): 
2.5  Hz  for  a  temperature  drop  of  30  K  and  for  the  very  high 
liquid  filling  1-a  =  0.95. 

DEFINITION  OF  PULSATING  DEVICES  FOR  TESTING 

A  first  step  in  a  logic  approach  to  develop  useful  pulsating 
two-phase  test  devices  is  to  dimension  these  such  that  test 
outcomes  can  directly  be  compared  to  the  presented 
performance  data  of  the  described  synchronised  oscillating 
single-phase  device. 

A  second  step  is  to  assume  that  the  single-phase  device 
(Fig.  11)  consists  of  85  identical  cylindrical  channels  with 
an  average  internal  diameter  d  =  1  mm.  This  represents 
the  actual  configuration  of  31  glass  capillaries  (1  mm  ID) 
and  54  triangular  channels  present  between  these 
capillaries,  thus  yielding  the  total  liquid  cross-sectional  area 
A  =  85*71/4  =  67  mm^.  For  simplicity  reasons  the  tube 
length  is  assumed  to  be  equal  to  the  displacement  length: 
L  =  Az  =  125  mm.  The  power  transported  by  each  of  these 
capillaries  can  be  calculated  from  the  data  presented.  For 
the  maximum  transport  case,  being  for  frequency  f  =  8  Hz 
and  temperature  difference  AT  =  56  K,  this  becomes  (n/4)* 
cf  *2.9*10®  =  (71/4)  *10’®»  2.3  W. 

The  third  step  is  the  simplification  to  consider  the  pulsating 
two-phase  device  (Fig.  14)  to  be,  in  essence,  a 
configuration  of  identical,  parallel  elements,  each  one 
transporting  the  same  amount  of  power,  driven  by  the 
vapour  pressure  difference  between  the  heat  input 
(evaporator)  section  and  the  cooiing  (condenser)  section. 
In  addition,  the  working  fluid  and  dimensioning  of  the  two- 
phase  and  single-phase  devices  are  identicai:  the  working 
fluid  is  water,  the  capiliary  diameter  d  =  1  mm,  the 
evaporator  and  condenser  length  are  Le  =  U  =  L  =  125 
mm.  In  the  first  approximation,  it  is  assumed  there  is  no 
adiabatic  section. 
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The  main  differences  between  the  two  devices  pertain  to 
the  driving  mechanism,  the  heat  transfer  processes  and 
the  heat  transfer  locations. 

A  mechanical  actuator  is  the  driver  of  the  osciiiating  axial 
movement  of  the  liquid  in  the  single-phase  device.  Only 
specific  heat  is  exchanged  over  the  entire  capillary  tube 
length  L  in  two  sequential  radial  (conduction)  steps.  In  the 
first  half  of  the  period,  heat  is  transferred  in  radial  direction 
from  the  hot  fluid  in  the  core  to  the  thin  Stokes  layer  (and 
the  tube  wall).  In  the  second  half  of  the  period,  this  heat  is 
moved  back  to  the  cold  fluid  brought  into  the  core. 

In  the  two-phase  device  heat  is  simultaneously  exchanged 
in  radial  direction,  mainly  by  conduction,  plus  some 
convection,  at  two  different  locations.  The  heat  is  fed  via 
the  wall  to  the  working  fluid  in  the  hot  input  (evaporator) 
section.  The  heat  is  extracted  from  the  fluid  via  the  wall  in 
the  cold  section  (condenser).  This  heat  transfer,  via  the 
specific  heat  of  the  liquid,  looks  more  or  less  identical  in  the 
two  systems.  The  transfer  difference  (i.e.  two-step 
sequential  at  one  location,  respectively  simultaneous  heat 
addition  and  extraction  at  two  different  locations)  suggests 
that  it  is  reasonable  to  assume  that  the  transported  power 
in  the  two-phase  case  is,  for  AT  =  56  K,  twice  the  value  for 
a  single-phase  capillary  (4.6  W).  As  an  alternative  it  can 
also  be  assumed  that  AT  =  28  K  only,  at  a  power  transport 
of  2.3  W. 

However,  there  is  an  additional  latent  heat  transfer 
contribution  in  the  two-phase  device:  The  heat  transported 
via  the  vapour  bubble  that  grows  in  the  heat  input  section 
(evaporator)  by  evaporation  of  a  liquid  micro-layer  (Ref. 
38).  This  bubble  collapses  in  the  condenser,  releasing  its 
latent  heat.  The  pressure  difference,  between  the  (super¬ 
heated)  vapour  in  the  evaporator  and  the  saturated  vapour 
in  the  condenser,  is  the  driving  force  moving  the  hot  liquid 
slug  from  evaporator  to  condenser,  plus  moving  at  the 
same  moment  a  similar  cold  slug  from  condenser  back  to 
the  same  or  a  neighbouring  evaporator.  The  power 
transported  by  latent  heat  can  be  obtained  by  calculating 
the  energy  needed  to  create  8  bubbles,  of  length  L  and 
diameter  d,  per  second.  Consequently  one  obtains 
8*(ji/4)*d^*L*pv*h|v  =  8*(7i/4)*10'®*  0.125  *  0.2  *  2.25*  10® 
=  0.45  W,  which  constitutes  a  non-negligible  contribution. 

The  pressure  head  across  the  capillary  single-phase  water 
system  can  be  calculated  as  follows.  The  displacement  of 
125  mm  at  frequency  8  Hz  yields  a  liquid  velocity  v  =  2  m/s. 
For  water  around  300  K,  the  Reynolds  number  Rei  is 
around  2000,  which  means  laminar  flow.  Consequently  the 
required  pressure  drop  is  8  kPa,  according  to  the  equation 

Ap  =  4*(16/Rei)(L/D)(py/2).  (39) 

In  the  corresponding  two-phase  device  the  required 
pressure  difference  has  to  be  far  larger,  because  of  several 
reasons.  In  the  first  place  twice  the  single-phase  device 
mass  (a  hot  and  a  cold  slug)  has  to  be  moved.  Secondly, 
this  double  mass  has  to  be  forced  through  a  180  degrees 
bend  instead  a  straight  channel.  Further,  the  length  of  the 
adiabatic  section  (La)  is  of  course,  in  reality,  never  equal  to 


zero.  Finally,  the  process  concerns  all  except  fully 
developed  flow,  hence  there  is  a  liquid  acceleration  term  to 
be  added. 

To  get  a  feeling  for  the  magnitude  of  these  pressure 
enlarging  effects,  the  length  La  is  taken  to  be  also  125  mm 
(La  =  L),  as  example.  This  has  impact  on  the  contribution 
of  the  power  transport  via  the  latent  heat  of  evaporation. 
This  contribution  will  be  around  0.9  W,  since  the  length  of 
each  of  the  8  vapour  bubbles,  being  generated  and 
collapsing  in  one  second,  is  La  -i-  L  (hence  2L,  instead  L). 
The  average  liquid  velocity  becomes  v  =  4  m/s.  For  water 
around  300  K,  the  Reynolds  number  Rei  now  lies  around 
4000,  which  means  turbulent  flow.  Consequently  the 
required  pressure  drop  has  to  be  calculated  according  to 

Ap  =  4  *  0.0791  Re’'^  (L / D)  (p,  v^/ 2).  (40) 

As  discussed  in  textbooks,  the  effect  of  the  two  bends  can 
be  accounted  for  by  adding  an  extra  length  of  50  D.  The 
pressure  drop  can  now  be  calculated,  according  to 
equation  (40),  by  inserting  the  different  parameter  values 
and  by  replacing  L  by  2*(2L  -i-  50  D).  The  result  is  480  kPa. 
To  be  complete  an  inertia  term,  accounting  for  the 
acceleration  of  the  slugs  eight  times  per  second,  has  to  be 
added:  8*(pi  v^/  2),  hence  16  kPa,  yielding  about  500  kPa 
for  the  pressure  difference  required.  Consequently,  it  can 
be  concluded  from  the  water  curve  in  figure  14,  that  the 
two-phase  device  has  to  operate  at  a  hot  section 
temperature  of  at  least  80  °C,  to  be  able  to  deliver  the 
pressure  drop  required.  The  figure  makes  also  clear  that 
more  or  less  comparable  powers  can  be  transported  by 
ammonia,  R12,  acetone,  etc.  Though  at  comparable  AT’s, 
this  will  be  realised  at  far  lower  operating  temperatures,  as 
these  fluids  show  a  steeper  dp/dT-relation. 

The  results  of  the  above  simple  approach  and  of  the 
detailed  modelling  of  the  physical  processes,  including 
mass-spring  simulations,  currently  is  compared  and  will  be 
compared  in  the  near  future  to  experimental  data,  resulting 
from  further  experimenting  at  NLR.  The  experimental 
activities  include  many  high-acceleration  experiments  on  a 
rotation  table.  The  experiments  have  been,  are,  and  will  be 
carried  out  both  with  all-metal  devices,  and  with  all-glass 
devices  (Ref.  36).  Additional  experiments  will  be  executed 
using  a  torus-like  or  spring-like  configuration  of  transparent 
(teflon  or  polyethylene)  flexible  tubing,  equipped  with  a 
simple  one-way  valve  to  influence  direction  and  frequency 
of  the  periodic  behaviour  of  a  closed-loop  configuration. 
Experiments  pertain  to  various  working  fluids  and  different 
locations  of  hot  and  cold  sections,  to  different  lengths  of 
adiabatic  section,  to  various  orientations,  and  to  different 
acceleration  levels  in  various  directions.  Results  between 
modelling  and  experimental  outcomes  will  be  subject  of  a 
publication,  which  is  being  prepared  for  presentation  early 
next  year  (Ref.  42) 

CONCLUDING  REMARKS 

In  conclusion  it  can  be  remarked  that  the  preceding  review 
summarises  the  approach  for  the  thermal-gravitational 
scaling  of  two-phase  systems.  Results  of  similarity 
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considerations  are  given.  Modeiiing  results  are  discussed, 
including  comparison  with  various  experimental  data.  The 
discussions  concern  not  only  gravity-assist  system  aspects, 
but  also  the  issues  of  operation  against  gravity  or  super¬ 
gravity.  The  systems  considered  include  capillary  pumped 
and  mechanically  pumped  two-phase  loops,  pulsating  two- 
phase  loops,  and  other  single-phase  and  two-phase 
oscillating  (pulsating)  devices.  A  simple  rationale  is 
described,  to  get  a  (quantitative)  feeling  to  be  used  for  the 
development  of  pulsating  devices  for  useful  experiments. 
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NOMENCLATURE 

A  area  (m^’ 

Ar  Archimedes  number  (-) 

B  proportionality  factor  (-) 

Bo  Bond  number  (-) 

Boil  boiling  number  (-) 

C  conductance  (W/K) 

Cp  specific  heat  at  constant  pressure  (J/kg.K) 
D  diameter  (m) 

d  diameter  of  capillary  or  of  curvature  (m) 

E  enhancement  factor  (-) 

E6  Eotvos  number  (-) 

Eu  Euler  number  (-) 

Fr  Froude  number  (-) 

f  frequency  (Hz) 

g  gravitational  acceleration  (m/s^) 

H  enthalpy  (J/kg) 

h  heat  transfer  coefficient  (W/m^.K) 

hiv  latent  heat  of  vaporisation  (J/kg) 

j  superficial  velocity  (m/s) 

L  length  (m) 

Ma  Mach  number  (-) 

Mo  Morton  number  (-) 


Mu  Inverse  viscosity  number  (-) 
rh  mass  flow  rate  (kg/s) 

Nu  Nusselt  number  (-) 

p  pressure  (Pa  =  N/m®) 

Pr  Prandtl  number  (-) 

Q  power  (W) 

q  heat  flux  (W/m^) 

Re  Reynolds  number  (-) 

T  temperature  (K  or  °C) 

r  radius  (m) 

S  slip  factor  (-) 

t  time  (s) 

V  velocity  (m/s) 

V*  dimensional  velocity  (-) 

We  Weber  number  (-) 

Wo  Womersley  number  (-) 

X  vapour  quality  (-) 

z  axial  or  vertical  co-ordinate  (m) 

a  vapour/void  fraction  (volumetric)  (-) 

p  constant  in  eq.  (13)  (-) 

5  surface  roughness  (m) 

A  difference,  drop  (-) 

K  thermal  diffusivity  (m^/s) 

A  thermal  conductivity  (W/m.K) 

p.  viscosity  (N.s/m^) 

v  angle  (with  respect  to  gravity)  (rad) 

7t  dimensionless  number  (-) 

p  density  (kg/m®) 

o  surface  tension  (N/m) 

Subscripts 

a  acceleration,  adiabatic,  axial 

c  condenser,  cold 

ch  (capillary)  channel 

e  evaporator 

eff  effective 

f  friction 

h  hot 

I  inner 

g  gravitation 

I  liquid 

m  momentum,  model 

0  reference  condition,  outer 

p  pore,  prototype 

r  radial 

s  entropy 

t  total 

tp  two-phase 

v  vapour 

w  water 


